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Abstract 

In  this  contribution  we  propose  the  use  of  general  bianisotropic  inhomogeneous  materials  for  both 
radiating  (i.e.  microstrip  antennas)  and  transmissive  (i.e  waveguides,  transmission  lines,  filters, 
etc.)  components  in  the  microwave  and  millimetre  wave  frequency  ranges.  The  analysis  of  such 
components  is  essentially  developed  via  a rigorous  full-wave  approach  in  the  spectral  or  in  the 
spatial  domain  depending  on  the  complexity  of  the  boundary  value  problem  to  be  solved  and  on 
the  numerical  technique  to  be  used.  For  the  transmissive  components  a fundamental  mode 
approximate  analysis  is  also  employed  in  order  to  obtain  very  fast  design  formulas. 


1.  Introduction 

Nowadays,  materials  with  a continuous  spatial  variation  of  the  refraction  index  can  be  easily 
synthesized  in  the  optical  field  for  graded  index  optical  fibers  and  planar  waveguides.  Also  in  the 
solide  - state  applications  we  can  control  the  carrier  concentrations  on  the  atomic  scale,  allowing,  thus, 
the  synthesis  of  components  with  an  essentially  continuous  conductivity  profile.  Since  in  the 
microwave  and  millimetre  wave  frequency  ranges  the  technological  processes  used  in  the  optical  and 
solid  - state  applications  cannot  be  effective,  the  synthesis  of  such  materials  could  be  realized  only  by 
means  of  artificial  dielectrics. 

Artificial  dielectrics  with  electromagnetic  properties  varying  with  the  spatial  coordinates  can 
be  synthesized  via  non-uniform  metallic  inclusions  in  a host  material.  In  addition,  by  using,  for 
instance,  omega  metallisations,  we  can  think  to  introduce  an  anisotropic  behaviour  together  with  a 
magneto-electric  effect.  In  such  a way,  thus,  we  can  obtain  a general  inhomogeneous  bianisotropic 
material. 

The  use  of  the  inhomogeneous  materials  has  been  recently  considered  in  some  works  [1-9]  and,  as 
shown,  leads  to  some  improvements  both  in  circuit  and  in  antennas  applications.  The  dielectric 
inhomogeneity,  for  instance,  can  be  used  to  control  the  filtering  and  matching  properties  of  a 
transmission  line  [1-4]  or  the  radiation  properties  of  microstrip  antennas  (directivity,  efficiency  and 
bandwidth)  [5-9]. 

The  analysis  of  such  components  has  been  developed  by  means  of  the  rigorous  full-wave  approach 
both  in  the  spectral  domain  associated  to  the  Method  of  Moments  (MoM)  and  in  the  spatial  domain 
associated  to  a variational  formulation  and  a Finite  Element  Boundary  Integral  (FE-BI)  numerical 
procedure.  Moreover,  in  certain  cases,  especially  when  the  propagation  characteristics  of  a component 
are  investigated  an  approximate  approach  based  on  the  TEM  or  quasi-TEM  propagation  can  be  used 
together  with  the  non-uniform  transmission  line  theory. 
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The  aim  of  this  contribution  is  to  extend  the  previous  works  [1-9]  to  the  case  of  inhomogeneous 
bianisotropic  materials  leading  to  a closed  form  representation  of  the  Green’s  dyads  in  the  spectral 
domain  and  to  a new  variational  formulation  in  the  spatial  domain. 


2.  Theory 

The  constitutive  relations  for  an  inhomogeneous  bianisotropic  material  are: 

[D  = e(r)  • E + £(r)  • H 
[B  = C(r)*E+|i(r)*H 


Assuming  a harmonic  time  dependence  exp[jcot]  the  Maxwell  equations  in  a source  free  region  can  be 
written  as: 


V x E = -jco^(r)  • E-  jco|i(r)  • H 

V x H = jcoe(r)  • E + jco^(r)  • H 

' V*[e(r)*E  + 4(r)*H]  = 0 (2) 

V*[C(r)»E+fA(r)*H]  = 0 

Since  the  analysis  of  radiating  and  transmissive  components  is  developed  in  different  ways  we  will 
present  it  in  two  separate  sub-sessions. 


2.1  Radiating  components 

Here  we  consider  planar  integrated  antennas  of  the  kind  depicted  in  Fig.  1 . 


(a)  (b) 


Fig.  1 Microstrip  patch  antenna  configuration:  Open  planar  integrated  patch  antenna  (a);  Cavity  backed  patch 
antenna  (b). 

The  electromagnetic  formulation  of  the  antenna  in  Fig. la  is  based  on  the  full-wave  analysis  in  the 
spectral  domain  and  the  numerical  solution  of  the  associated  boundary-value  problem  is  obtained  via 
the  method  of  moments  (MoM).  The  fundamental  step  in  this  approach  is  the  derivation  of  the  spectral 
Green’s  dyad  for  the  planar  structure.  In  this  case  we  consider  the  variation  of  the  electromagnetic 
parameters  only  along  the  stratification  direction  and,  starting  from  the  curl  Maxwell  equations,  we 
obtain,  after  a change  of  the  reference  system  [6],  the  following  relations  in  the  spectral  domain: 
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^^  = Cvv(y).V(y)+CVI(y).I(y) 
dy 

® = Cjv  (y)  • V(y) + Cn  (y)  • I(y) 


where  V = [Ev,  Eu]T,  I = [HU,  Hv]Tand  the  elements  of  the  matrices  CVY(y),CVI(y),CIV(y),C1I(y) 

are  complicate  ratios  involving  the  electromagnetic  parameters.  In  order  to  decouple  system  (3)  we 
impose  Cvv(y)  = 0 and  Cn(y)  = 0 obtaining,  thus,  after  three  derivative  operations,  the  following 
fourth  order  decoupled  differential  equations: 


d4Ev 

dy4 

d4Eu 

dy4 


+ Av3(y)-^  + Av2(y)-^  + Avl(y)^+Av0(y)Ev=0 

d3F  d2E  dE 

+AU3(y)^  + Au2(y)^  + Aul(y)^-+Au0(y)Eu=0 


where  the  non  constant  coefficient  Au(y)  and  Av(y)  depend  on  the  elements  of  the  matrices  Cjj(y). 
Once  the  solutions  of  such  equations  are  found,  the  overall  spectral  electromagnetic  field  is  also 
known  in  the  grounded  inhomogeneous  slab  and  the  Green’s  dyad  can  be  easily  determined. 

The  application  of  the  MoM  in  order  to  calculate  the  unknown  current  density  on  the  patch,  is  now  a 
straightforward  matter  [7-8], 

On  the  other  hand  the  electromagnetic  formulation  for  the  antenna  in  Fig.  lb  is  based  on  a 
variational  approach  together  with  the  implementation  of  the  Finite  Element  Boundary  Integral  (FE- 
BI)  numerical  method.  The  functional  F(E)  that  has  to  be  stationary  for  arbitrary  variations  of  the  E 
field  is  given  by: 


F(E)  =<  jx(r)a  1 »(VxEa  +jco|(r)«Ea),VxE  + ja)^(r)»E>+a)2<Ea,e(r)»E>- 

-<Ea,j©J>-< j(0ja,E>+2(0—  f(nxEa)*»  f G(r,r>rzxE(r')ldS’ 

"Ho  si  L ap 


dS 


(3) 


For  the  meaning  of  the  apex  a,  of  the  vector  J,  of  the  surface  Sap,  please,  refer  to  [5]. 

After  discretizing  the  cavity  into  finite  elements  of  proper  shapes  (in  this  case  bricks  are  used), 
and  after  expanding  the  electric  field  by  using  the  edge-based  vector  basis  functions  as  in  [5]  we 
obtain  an  algebraic  system  to  be  solved  of  the  form  A • X = B . The  solution  of  such  a linear  system 
returns  the  electromagnetic  field  in  the  cavity  and  on  the  aperture  surface  and,  thus,  the  derivation  of 
the  antenna  characteristics  is,  now,  a straightforward  matter. 


2.2  Transmissive  components 

Here  we  consider  both  planar  and  not  planar  waveguides  as  shown  in  Fig.2.  These  waveguides  have 
been  studied  by  means  of  different  approaches.  When  the  variation  of  the  electromagnetic  parameters 
is  only  along  the  stratification  axis  striplines  and  microstriplines  have  been  analyzed  in  the  spectral 
domain  by  using  the  same  approach  previously  shown  simply  introducing  the  z dependence  exp[-jkzz] 
in  the  Maxwell  equations.  The  integral  equation  written  on  the  strip  can  be  solved  again  via  the  MoM 
leading  to  an  eigenvalue  algebraic  system  whose  solution  returns  the  propagation  characteristics  of  the 
waveguides.  Striplines,  microstriplines  and  rectangular  waveguides  with  arbitrary  variations  of  the 
parameters  only  on  the  cross  sections  have  been  analyzed  by  means  of  the  full-wave  approach  in  the 
spatial  domain  together  with  the  Method  of  Lines  (MoL)  numerical  procedure  [10],  Finally,  when  the 
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parameter  variation  is  only  along  the  power  flow  direction,  for  the  striplines  and  microstriplines  the 
approximation  of  a TEM  or  quasi-TEM  propagation  has  been  used  and  the  solution  has  been  obtained 
by  means  of  the  non  uniform  transmission  line  theory  as  in  [1-4]. 


(a)  (b)  (c) 

Fig.  2 Three  different  kinds  of  waveguides  here  considered:  Microstripline  (a);  Stripline  (b);  Rectangular 
Waveguide  (c). 


3.  Conclusion 

In  this  contribution  we  have  proposed  the  use  of  general  bianisotropic  inhomogeneous  materials  for 
both  radiating  and  transmissive  components  in  the  microwave  and  millimetre  wave  frequency  ranges. 
The  analysis  of  such  components  has  been  essentially  developed  via  a rigorous  full-wave  approach  in 
the  spectral  or  in  the  spatial  domain  together  with  different  numerical  techniques  (MoM,  FE-BI,  MoL, 
etc.).  Particularly,  a new  variational  formulation  for  cavity  backed  microstrip  patch  antennas  loaded 
with  general  inhomogeneous  bianisotropic  dielectrics  and  new  analytical  expressions  of  the  spectral 
dyadic  Green’s  functions  for  grounded  slab  loaded  with  the  same  materials  have  been  derived  as 
extensions  of  works  previously  published. 
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